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£f) ' corrections in QCD and field theory. 
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' 1. Introduction 

t-H ! 

The precision of the experimental data on electroweak interactions and 
QCD is now very high and it is expected to become significantly higher 
| within the next few years. This has triggered a substantial refinement in 

the corresponding theoretical calculations. Yet, already now for certain 
experimental quantities the theoretical uncertainty is one of the major open 
' questions in the interpretation of the data and in the search for signals of 

C^-i, physics beyond the Standard Model. A striking example is the need for 

O j' a precise determination of the gauge couplings at the weak scale, which is 

the prerequisite for investigation of possible unification of couplings at some 
GUT scale. 

One of the reasons for this current state of affairs in the relation be- 
tween the theory and experiment is that computation of high orders in 
' perturbation theory for quantum field theories, and especially non-abelian 

gauge theories in 3+1 dimensions is extremely hard. State-of-the-art calcu- 
lations available today for this kind of theories have reached, after a very 
large effort, the 3-rd and the 4-th order in a s , for observables and for the 
/3-function, respectively [jl|, Q ||. Without a major breakthrough in the 
relevant techniques it is unlikely that exact results for the next order will 
become available in the foreseeable future. Moreover, even if explicit ex- 
pressions for very high order terms do become available, we still have to 
deal with the fact that the perturbative series of interest are asymptotic, 
with zero radius of convergence and usually are not even Borel summable. 
In this talk I will review an approach which has been recently suggested to 
deal with some of these problems. 
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2. Perturbation Theory: Diseases and a Promising Therapy 

As mentioned in the Introduction, the perturbation series in QCD is 
expected to be asymptotic with rapidly growing coefficients: 

oo 

S(x) = Y c n x n , x = — , c n ~ n\KV (1) 

for some coefficients K, 7 0,0]. Anyone who wants to make use of QCD 
perturbation theory to carry out precision analysis of observables has to 
face several practical problems: 

• only few first orders in ([l]) are known for any observable (n < 3) 

• the series has zero radius of convergence 

• the series is usually not Borel summable. Borel summation is a trick 
that sometimes works for summing series with factorial divergence. 
Consider the series for S(x) in eq. (|l]). We can define a new function, 
S(y), whose series is obtained from (|l]) by dividing the n-th term by 



n!, 



00 



,1 

n=0 



<%) = £ 3 h n (2) 



If the new series is convergent, the original function S(x) can be ob- 
tained by the so-called inverse Borel transform, 

S(x) = - / e-y' x S(y)dy (3) 



x Jo 

provided <S(y) has no singularities along the integration path. 
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Unfortunately, in QCD it is 
known that for a generic ob- 
servable S(y) has poles on 
both the positive and nega- 
tive y axis in the complex y 
plane. These are usually re- 
ferred to as infrared and ul- 
traviolet renormalons, respec- 
tively. The figure to the left 
shows a schematic description 
of poles in the Borel transform 
of a generic series for a QCD 
observable. 
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The presence of singularities along the integration path makes the 
integral (||) ill-defined. One can try to define it by going around 
the poles, but this introduces an ambiguity proportional to the pole 
residue, since different deformations of the integration path will give 
different results. 

• renormalization scale dependence: finite-order perturbative predic- 
tions depend on the arbitrary renormalization scale /i through the 
coupling, a s = a s (/x). This renormalization scale is most pronounced 
at leading order in perturbation theory and decreases with the inclu- 
sion of higher order terms. 

• renormalization scheme dependence: in principle, the theory can be 
renormalized in any valid renormalization scheme, yielding the same 
predictions for any physical observable. In practice, when we work 
with a finite number of perturbative terms, the results depend on the 
renormalization scheme. 



There is no "miracle cure" which would solve these problems completely. 
However, we should and can minimize their effect. Thus the practical issue 
is how to get the best possible precision, given a fixed number of terms in the 
perturbative expansion. In the following, I will discuss one method which 
has already shown considerable progress towards this goal. The method is 
based on the so-called Pade Approximants (PA-s) 0-[17]. 



3. Pade to the Rescue 



Pade Approximants |jj8| , 19 1 are rational functions chosen to have a Tay- 
lor expansion equal the perturbative series to the order calculated. Given a 
series 

S(x) = c + c lX + c 2 x 2 + c 3 x 3 + • • • + c N x N + 0{x N+l ) (4) 
one can always find a rational function 

such that [L/M] has the Taylor series 

[L/M] = c + c x x + c 2 x 2 + c 3 x 3 + ■■■ + c N x N + c N+1 x N+1 + ■■■ (6) 
The rational function in (||) is called the [L/M] Pade Approximant. 
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It is important to keep in mind that at a given finite order in x the 
[L/M] Pade in (|5|) is formally as valid representation of S(x) as the original 
perturbation expansion. Moreover, in practice the PA-s turn out to posses 
many important and useful properties which are absent in the straightfor- 
ward perturbation theory. 

Thus, even though PA is constructed to reproduce the series (|j) only up 
to order N, it turns out that under rather mild conditions the next term in 
the Taylor expansion of the PA in eq. (||), cn+i, provides a good estimate, 
c^+i = cjv+i. We call it the Pade Approximant Prediction (PAP), of the 
next coefficient cn+i in the series (0): 



est 
C N+1 _ , 



« 1 (7) 



CN+1 

and for sufficiently large N the relative error decays exponentially fast, 



est 
N+l 



1 



CN+1 



~ e~ aN , a = const (8) 



Let us consider some simple examples, starting with the trivial case of a 
single-pole geometric series 

A 00 

^ = <») 

n=0 

It is easy to convince oneself that in this case the [L/M] Pade is exact for 
L > 0, M > 1. For example, if we attempt to construct a [10/10] Pade of 
(H), we will find that the a priori 10-th degree polynomials in numerator 
and denominator reduce to a degenerate case of a constant and 1-st degree 
polynomial, respectively, 

10/10 = = io 

Qio{x) 1 - Bx 

Once this is clear, the extension to a sum of finite number of poles in obvious, 

K A oo 

; — = Cn%n =^ [ L / M ] exact for L > K - 1, M > K 

i=i 1 ~ BiX n=0 

One can also show that for an infinite number of isolated poles, i.e. when 

oo 

f(x) = c n x n is a meromorphic function, the sequence of [L/L + k] for k 
o 

fixed converges to f(x) as L — > oo, 

[L/L + k] > f(x); k = 0,±1,±2, ... (12) 
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A somewhat less intuitive, but very important result is that in certain cases 
the Pade sequence [L/L+k] converges exponentially fast in L to the correct 
function even for a factorially divergent asymptotic series with zero radius 
of convergence. A classical example [19] is the function 



1 + xt 



nl 



(13) 



Here again it turns out that [L/L + k] — ► g(x) as L — > oo. 

The crucial property of the series in ( Jl3| ) which makes this possible is 
that it has alternating signs. It is easy to show that this implies that all 



the poles of the Borel transform of (13) are on the negative real axis, and 
hence that the series is Borel summable. More generally, when the series is 
Borel summable, Pade will converge to the correct result. 

It is interesting to note that the exponentially fast convergence of PA-s 
is not limited to meromorphic functions. As a simple example, consider 
the hypergeometric function F(— 1 |, l,x) [2C], which has a cut for x > 1. 
Fig. H shows that despite the cut, the diagonal PA-s evaluated at x = 0.2 
converge exponentially fast. 



Pade approximants to F(- 1/2,3/2, 1 ,x), x = 0.2 
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Fig. 1. Exponentially fast convergence of Pade despite the presence of a cut: rela- 
tive error of diagonal [L/L] Pade for the hypergeometric function F(— 1/2, 3/2,1, x) 
which has a cut for x > 1. The Pade is evaluated to the left of the cut, at x = 0.2. 
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4. Applications to Quantum Field Theory 



While such mathematical examples are instructive, in order to gain con- 
fidence in the method, we need to see how it fares on high-order series taken 
from quantum field theory. As the first test case, we consider the scalar field 
theory with Gaussian propagators. High-order perturbation expansions of 
Green's functions in this theory have been computed in Ref. [21|. Fig. [2| 
demonstrates the convergence of PAP for the relevant coefficients for the 
4-point Green's function in D = 4 |j. The relative error is ~ 10" 3 at 5-th 
order. For comparison also shown are relative errors of estimates based on 
asymptotic behavior of large orders in perturbation theory, as given in || . 
Clearly, at 5-th order Pade does better by about 4 orders of magnitude. 



Predicting the next term in P.T. for \(p in D=4: 
Pade vs. large — order asymptotic formulae 
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Cj xact are tor 4-Point Green's Function, Bervillier et al., PRD, 17,2144(1978) 

Fig. 2. Pade predictions (crosses) for high-order terms in the perturbative expan- 
sion for 4-point Green's function in a scalar field theory in D = 4. For comparison 
also included are predictions based on asymptotic large-order behavior (circles). 

If information is available about the asymptotic behavior of c n , it is 
possible to obtain an explicit expression for the the error formula on the 
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r.h.s. of (||). For example, we have demonstrated that if 

ee - 1 * ±, (14) 

C n+1 n 

as is the case for any series dominated by a finite number of renormalon 
singularities, then <5[l/m] defined by 

8{L/M] = (15) 

CL+M+1 

has the following asymptotic behaviour 

Ml 

S [l/m]^-^m, where K = L + M + aM (16) 

and where a is a number of order 1 that depends on the series under con- 
sideration. For large L, M eq. (|l6|) yields an exponential decrease of the the 
error, as in eq. (||]). 

This prediction agrees very well with the known errors in the PAP's § 
for the QCD vacuum polarization D function calculated in the large N* 
approximation [p^j , as seen in Fig. |3|a. 

One can repeat this exercise also for the Borel transform of the D func- 
tion series. As mentioned earlier, a generic Borel transform is characterized 
by the presence of poles (more generally, branch points) on the real axis. 
In view of this, we expect an even faster convergence in this case, since the 
Pade, being a rational function, is particularly well-suited to reproduce this 
analytic structure. Indeed, it turns out that in this case e n ~ 1/n 2 and 

(Ml) 2 

${M/M] ^ ~ -J^M (17) 

which agrees very well with the corresponding PAP results shown in Fig. [3|b 
1,0. 

The high degree of agreement between the analytical error estimates 
in eqs. ( |I6| ) and ( |i~7| ) and the actual errors in PAP suggest that one can 
substantially improve the PAP method by systematically including the error 
estimates 6<l /m] cls a correction, yielding the Asymptotic Pade Approximant 
Predictions (APAPs): 

^APAP _ C L+Af+1 ri oA 
CL+M+1 ~ 1 + 5L+M+1 (18) 



where cf^ M+1 is the original PAP prediction without the additional correc- 



tion, as in e q. (|i~5| ), and <5l+m+i is obtained by fitting (16) to the known 
lower orders flU 
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Fig. 3. Relative errors in the [L/M] Pade Approximant Predictions |7|, O] (a) for 
the QCD vacuum polarization D-function series, evaluated to all orders inthc large- 
Nf approximation [ p2[ (the rate of convergence agrees with expectations for a series 
with a discrete set of Borel poles), and (b) for the Borel transform of the D-function 
series, where the convergence is particularly striking. The straight lines correspond 
to the exponential decay given by the respective error formulae, eqs. ([l6]) and (|l7|). 
The crosses, diamonds and squares correspond to M = L + l, L + X,L, respectively. 



The APAP method results not only in a substantial improvement of 
the PAP estimates, but also significantly reduces the difference between the 
predictions based on different [L/M] values at a given order, L + M=fixed. 

The method has been applied to the Bjorken sum rule for the difference 
of first moments of proton and neutron structure functions gi(x,Q 2 ) in 
polarized deep inelastic scattering [23]. For Nf = 3 the sum rule reads 



\g p 1 {x,Q 2 )-g n l {x,Q 2 )] = l -\g A \ 
o o 



1 — x — 3.58x 2 - 20.22x 3 + c 4 x 4 + . . . 



where x = a(Q 2 )/ir and the exact expression for C4 is still unknown. 



(19) 
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The PAP and the corresponding APAP estimates of C4 are 
[2/1]: c PAP ^-lU — » ct PAP ^-131 

[1/2]: c£ AP «-lll — » cf PAP «-130 



(20) 



Clearly, APAP estimates show significantly less spread than the correspond- 
ing PAP estimates. Remarkably, the APAP estimates of £4 show an almost 
perfect agreement with an independent estimate, based on a completely 
different method @: c 4 = -130 ! 

As already mentioned, a typical finite order perturbative series such 



as (19) exhibits a spurious renormalization scale and scheme dependence. 



Schematically we have, 

Six) = en + c\x + C2X 2 + C3X 3 + C4X 4 + C5X 5 + • • • (21) 
v , ' , ' 

physical observable: exactly known partial sum, unknown 
scale and scheme scale and scheme dependent higher order terms 

independent 

Replacing a finite-order perturbative series by a Pade is equivalent to adding 
an infinite series of estimated terms generated by the rational approximant. 
If such an estimate is accurate, we expect to see a reduction in the renormal- 
ization scheme and scale dependence. As shown in Fig. ^, this expectation 
is fully realized when Pade is applied fl2[| to the Bjorken sum rule series in 
eq. ©. 

It turns out that this dramatic reduction in the scale and scheme depen- 



dence can also be understood on a deeper level. In Ref. [jl3[ it was shown 
that in the large- (3q limit, i.e. when the (5 function is dominated by the 
one loop contribution, the scale dependence is removed completely. This is 
because in this limit the renormalization scale transformation of a s reduces 
to a homographic transformation of the Pade argument. Diagonal PA's are 



invariant under such transformations [18]. Non-diagonal PA's are not to- 
tally invariant, but they reduce the RS dependence significantly [pi. In 
the real world the usual MS function includes higher-order terms beyond 
/3o- Still, in QCD with 3 < Nj < 5, the 1-loop running of the coupling 
is dominant and therefore PA's are still almost invariant under change of 
renormalization scale. 

A further related interesting development is the observation |l6| that 
the Pade approximant approach for resummation of perturbative series in 
QCD provides a systematic method for approximating the flow of momen- 
tum in Feynman diagrams. In the large-/?o limit, diagonal PA's generalize 



the Brodsky-Lepage-Mackenzie (BLM) scale-setting method [25] to higher 



orders in a renormalization scale- and scheme-invariant manner. 
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Bjorken effective-charge - a NNLO partial-sum 




Bjorken effective-charge - [0/2] PS 




Fig. 4. Bjorken effective charge for Q 2 = 20GeV 2 plotted as a function of the 
renormalization scale and scheme, as parametrized by the coupling x = a s {ii 2 )/ir 
and the second coefficient of the j3 function: C2 = /?2/A)- The upper plot shows the 
NNLO partial sum, while the lower plot shows the [0/2] Pade approximant. 
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5. Predicting the QCD (3 function at 4 and 5 loops 

Although no QCD observables have been calculated exactly beyond 
0(a 3 ), in fall of 1996 we had learned that a calculation of the 4- loop con- 
tribution to the QCD (3 function was under way, and likely to be published 



soon. The prediction of the unknown 4-loop coefficient [14| was therefore 
an important challenge and excellent testing ground for the new APAP 
method. 

As a warm-up exercise one can test the APAP method on the 4-loop (3 
function of the $ 4 theory with O(N) global vector symmetry p6|] , the latter 
being analogous to the SU(Nf) global symmetry of QCD. The results [14] 
for (3% in O(N) <J> 4 theory are shown in Figure ||. Clearly, the variant of 
APAP method denoted as {A)/n (see (14] for details) is markedly superior 
to the naive PAP. The 5-loop (3 function in this theory is also known [27, 28] 
and the corresponding APAP estimates also turn out to be very precise [|14| . 
Consequently this was the method of choice for the QCD 4-loop (3 function. 
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Fig. 5. The 4-loop /3-function coefficient (3-$ in $ 4 theory with O(N) symmetry. 
The exact results are denoted by black dots, joined by a solid line to guide the 
eye. Naive PAP results are denoted by diamonds, and APAP results obtained from 
the {A) jn type of correction are denoted by crosses. For comparison, also shown 
are APAP results obtained from the (B)/n 2 type of correction, denoted by open 
circles. 
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The strategy for computing (3 3 , the 4-loop (3 function coefficient, is as 
follows. We recall that (3 3 is a cubic polynomial in the number of flavors 

p 3 = A 3 + B 3 N f + C 3 N 2 f + D 3 N 3 f , (22) 

where D 3 = 1.49931 (For Nc = 3) is known from large- Nf calculations p9| . 
The known exact expressions for the 1-, 2- and 3- loop (3 function are used 
as input to APAP, to predict the value of (3 3 for a range of Nf values. The 
predictions for A 3 ,B 3 , C 3 are then obtained from fitting the APAP results 
for < Nf < 4 to a polynomial of the form (|22|) . 

Shortly after the APAP prediction of (3 3 |14| the exact result was pub- 
lished in Ref. One important lesson from the exact results is that they 
contain qualitatively new color factors, corresponding to quartic Casimirs, 
analogous to light-by- light scattering diagrams in QED. Such terms are not 
present at 1-, 2- and 3-loop level, and therefore cannot be estimated using 
the Pade method. Numerically the correction due to these new color factors 
is not very large, but in principle the PAP estimates should be compared 
with the rest of the exact expression, as is done in the first three columns 
of Table I. 
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Table 1. Exact four-loop results for the QCD (3 function, compared with the orig- 
inal APAP's in the first column, and improved APAP's obtained from a weighted 
average over negative Nf (WAPAP), as discussed in fj~7| . The numbers in paren- 
thesis are the error estimates from [[l4] . 

The APAP estimates for ^3, B 3 and C 3 seem quite satisfactory, until 
one realizes that the A 3 and B 3 terms in (|22|) have opposite signs and 
their magnitude is such that they almost cancel each other at Nf ~ 4. This 
means that for numerical prediction of (3 3 as function of Nf in the physically 
interesting range < Nf < 5 a better precision is required. Fortunately, 
it is possible to obtain such precision. This is accomplished by formally 
using negative values of Nf in the fitting procedure, so that no cancellation 
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occurs, and making a careful choice of the range of negative values of Nf 
used for the fit. Once the values of A3, B3 and C3 are obtained this way, 



one can use them to compute the physical predictions at positive Nf []17 1 . 
This procedure has been referred to as WAPAP, for "weighted APAP" . The 
corresponding results are shown and compared with exact results in the last 
two columns of Table I. We see a dramatic improvement in the precision. 

Figure || displays graphically predictions for /?3, as a function of Nf 
for the most interesting case Nc = 3. We plot the percentage relative er- 



rors obtained using various APAP-based estimation schemes [14, 17]: naive 
APAP's fitted with positive Nf < 4 (diamonds), naive APAP's fitted with 
negative Nf > —4, WAPAP's compared to the exact value of (3% includ- 
ing quartic Casimir terms, and WAPAP's compared to without quartic 
Casimir terms (crosses). We see that the latter are the most accurate for 
/? 3 in QCD. 
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Fig. 6. Predictions for ^3, as function of Nf, for Nc = 3. The percentage rel- 
ative errors are obtained using various APAP-based estimation schemes: naive 
APAP's fitted with positive Nf < 4 (diamonds), naive APAP's fitted with negative 
Nf > —4, WAPAP's compared to the exact value of 03 including quartic Casimir 
terms, and WAPAP's compared to /?3 without quartic Casimir terms (crosses). 
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In Figure [7| we show the error in the WAPAP prediction for (3% as a 
function of Nf, for Nc =3, 4, 5, 6, 7 and 10, once again omitting quartic 
Casimir terms from the exact result. The accuracy of these predictions is 
our best evidence for believing in the utility of the WAPAP method. 



1.00 
0.10 
0.01 

1.00 
0.10 
0.01 

1.00 
0.10 
0.01 

1.00 
0.10 
0.01 

1.00 
0.10 
0.01 

1.00 
0.10 
0.01 



2 Ntr 3 



Fig. 7. The percentage relative errors in the WAPAP prediction for 03 (compared 
to the exact result with quartic Casimir terms omitted), plotted vs. Nf for Nc =3, 
4, 5, 6, 7, 10. 



The WAPAP method does very well on the four-loop QCD (3 function, 
but the details of the method were fine-tuned after the exact results became 
available. In Ref. [|17]] predictions were also given for yet unknown 5-loop 
j3 function in QCD and 4- and 5-loop (3 function in N = 1 supersymmetric 
QCD. It is extremely interesting to see how well our predictions will do for 
these quantities. 
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